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We present a study of the effects of temperature upon the excitation frequencies of a Bose-Einstein 
condensate formed within a dilute gas with a weak attractive effective interaction between the atoms. 
We use the self-consistent Hartree-Fock Bogoliubov treatment within the Popov approximation and 
compare our results to previous zero temperature and Hartree-Fock calculations The metastability 
of the condensate is monitored by means of the I — excitation frequency. As the number of 
atoms in the condensate is increased, with T held constant, this frequency goes to zero, signalling 
a phase transition to a dense collapsed state. The critical number for collapse is found to decrease 
as a function of temperature, the rate of decrease being greater than that obtained in previous 
Hartree-Fock calculations. 

PACS Numbers: 03.75.Fi, 05.30.Jp, 67.40. Db 

Mean field theories of the Bose-Einstein condensation of trapped alkali vapours have been extremely successful both 
qualitatively and quantitatively in determining the excitation frequencies of the condensates, especially at relatively 
low temperatures (< 0.7 Tc) pi. These calculations have been based upon the Popov approximation to the Hartree- 
Fock Bogoliubov (HFB) treatment, where the anomalous average of the fluctuating field operator is neglected [||. 
In all cases the study has been of alkali vapours with positive s-wave scattering lengths (i.e., repulsive effective 
interactions). The case of attractive interactions (^Li for example, as used in experiments at Rice University M) has 
not been treated in this manner. Calculations have, rather, been based on the zero temperature Gross-Pitaevskii 
equation (GPE) M or upon a Hartree-Fock variational calculation M. 

There are two main reasons why the HFB formalism was not used in the Hartree-Fock study referred to above. 
Firstly, in the case of negative scattering length, the HFB-Popov collective excitations of a homogeneous system are 
unstable at long wavelengths. Houbiers and Stoof [|| therefore found it more appealing to use the Hartree-Fock 
method, which has stable excitations at long wavelengths. In the case of the trapped gas this does not present a 
problem, as one is saved from the infra-red limit by the finite zero-point energy of the trap. From an alternative 
viewpoint, the finite size of the condensate eliminates very long wavelength excitations. The HFB-Popov theory is 
hence quite applicable for trapped gases. 

Secondly, there is the possibility that atoms with an attractive effective interaction can undergo a BCS-like pairing 
transition Q. This possibility is in fact included in the full theory that Houbiers and Stoof develop. However, in their 
numerical calculations , they ignore the possibility of pairing and the results presented are based on a Hartree-Fock 
treatment of Bose-Einstein condensation alone. If one is going to assume that their is no BCS transition, then a better 
description would appear to be that of the HFB-Popov formalism. This is the treatment adopted in this letter. 

The purpose of the present investigation is to determine the stability of the condensate against mechanical collapse, 
and the effects thereon of thermal excitations. It has been shown [Q in the homogeneous limit that the condensate 
is unstable at the densities required for BEC. In the trap, the additional kinetic energy can stabilise the condensate 
and a metastable state is possible. This state decays on a timescale which is long compared to the lifetime of the 
experiment, but only exists for condensates below a certain size. At some critical condensate number the condensate 
becomes unstable and collapses. This instability is characterised by the monopolar collective excitation going soft Q] 
(viz., the excitation frequency goes to zero). Various predictions for the critical number, Nc, have been made at zero 
temperature using the GPE and at finite temperatures using the Hartree-Fock treatment. Here we investigate the 
effects of temperature upon the collapse via the HFB-Popov approach as described briefly below. 

We make the usual decomposition of the Bose field operator into condensate and noncondensate parts; ip{r) = 
^(r) -I- '0(r). The condensate wavefunction (f>(r) is then defined within the Popov approximation by the generalised 
Gross-Pitaevskii equation (GPE) 
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Here, no(r) = |$(r)p and n(r) = {ip'' {r)ip{r)) are the condensate and noncondensate densities respectively. The Popov 
approximation p-no| consists of omitting the anomalous correlation {ip{r)^{r)), but keeping h{r). The condensate 
wavefunction in EqT^) is normalised to Nq, the total number of particles in the condensate. T4a;t(r) is the external 
confining potential and g = Anh a/m is the interaction strength determined by the s-wave scattering length a. For 
'^Li the value of a used is -27.3 Bohr radii. The condensate eigenvalue is given by the chemical potential ^ Q. 

The usual Bogohubov transformation, tplr) = ^i[uii'r)ai — i;*(r)dj], to the new Bose operators di and a] leads to 
the coupled HFB-Popov equations H 

£ui{r) - gno{r)vi{r) = EiUi{r) 

Cvi{v) - gno{r)ui{r) = -E,Vi{Y) , (2) 

with £ = — V^/2TO + Vea;t(r) + 2(77i(r) — /i = ft,o+.9'^o(r)- These equations define the quasiparticle excitation energies Ei 
and the quasiparticle amplitudes Ui and Vi. Once these quantities have been determined, the noncondensate density 
is obtained from the expression |q| 

i 

= ni(r) +n2(r), (3) 

where ni(r) is that part of the density which reduces to the quantum depletion of the condensate as T ^ 0. The 
component ri2(r) depends upon the Bose distribution, No{E) = {e^^ — 1)^^, and vanishes in the T ^ limit. 

Rather than solving the coupled equations in Eq. (y) directly, we introduce the auxiliary functions ipi (r) = 
Ui(r) ± Vi(r) which are solutions of a pair of uncoupled equations (a more detailed discussion of the method is 
presented in Hutchinson et al. of Ref.[l]). The two functions are related to each other by /loV'i — Eiijj^'. We 
note that the collective modes of the condensate can be shown to have an associated density fluctuation given by 
5ni(v) oc ^{r)'il}\~\Y). 

To solve these equations we introduce the normalised eigenfunction basis defined as the solutions of hQ(f)a{r) = 
£a4'a{^) a-iid diagonalize the resulting matrix problem. The lowest energy solution gives the condensate wavefunction 
<3?(r) = v^/Vo0o(r) with eigenvalue Eq — Q. 

The calculational procedure can be summarised for an arbitrary confining potential as follows: Eq.(|l|) is first solved 
self-consistently for $(r), with n(r) set to zero. Once $(r) is known, the eigenfunctions of Kq required in the expansion 
of the excited state amplitudes are generated numerically. The matrix problem is then set up to obtain the eigenvalues 
£■,;, and the corresponding eigenvectors Ca are used to evaluate the noncondensate density. This result is inserted 
into Eq.(y) and the process is iterated, keeping the condensate number iVo and temperature T fixed. The level of 
convergence is monitored by means of the noncondensate number, N and the iterations are terminated once N is 
within one part in 10^ of its value on the previous iteration |13]. In this way, we generate the self-consistent densities, 
no and h, as a, function of Na and T. 

We consider first the case of T = for an isotropic harmonic trap with a frequency equal to the geometric average 
of the frequencies corresponding to the Rice trap H, i^ = 144.6 Hz. This is the geometry considered previously by 
Houbiers and Stoof g and with whom we find qualitatively agreement. There are several signatures of a collapse 
of the condensate with increasing condensate number Nq. First, we can look at the behaviour of the convergence 
parameter (the total number of particles in the noncondensate for a given condensate number and temperature) used 
to monitor the convergence of the solution to the HFB-Popov equations. In Fig. 1 we show A^ as a function of iteration 
number for the three values A^o — 1243, 1244, and 1245. The convergence is clear in the first two cases, whereas in 
the final case the algorithm diverges catastrophically and no stable solution can be found. We therefore identify the 
critical number, Nc, of atoms in the condensate as 1244, beyond which the condensate is no longer metastable, but 
unstable to the formation of a dense solid phase. This value of Nc is slightly greater than the value of 1241 obtain by 
Houbiers and Stoof using the Hartree-Fock approximation. A second, more physical indicator of the collapse is the 
observed strong dependence of the excitation frequencies on the number of condensate atoms. In particular, we find 
that the I = mode goes soft as A^o approaches the critical number found above. We shall focus on this criterion for 
the instability in the following. 

We next consider a trap with confining frequency 150 Hz. The excitation frequencies are again calculated as a 
function of the number of particles in the condensate, both at T = and at finite temperature. The lowest lying 
modes at temperatures of 0, 200, and 400 nK are shown in Fig. 2. The lowest mode is the I — 1 Kohn mode, 
which corresponds to a rigid centre of mass motion. For a harmonic trap the excitation frequency of this mode 
should be identically equal to the trap frequency. However, the dynamics of the noncondensate are neglected in this 



treatment and the calculated excitations are those of the condensate alone, moving in the effective static potential 
Veff = Vext + 2gh{r). Due to the presence of the noncondensate, the effective potential is not parabolic and hence 
the generalised Kohn theorem does not apply. The Kohn theorem is approximately obeyed for low temperatures and 
low particle numbers since the noncondensate is either small, or relatively uniform over the extent of the condensate 
and hence docs not introduce a significant anharmonicity. It is only for higher temperatures near Nc where the 
noncondensate density is both large and sharply peaked around the centre of the trap that there is a marked deviation 
from the trap frequency. 

As mentioned above, the softening of the I — breathing mode is a signature of the instability from a metastable 
condensate to a completely collapsed state. For T = the critical number, N^ is found to be 1227, which is slightly 
lower than that obtained in the previous case with a stronger confining potential. This is the change in the critical 
number expected B on the basis of the dependence N^ oc 1/y^u^, which shows that the critical number increases as 
the trap confinement is relaxed. With increasing temperature, the frequency of the I — mode is found to go to zero 
at lower condensate numbers. This is because the attractive nature of the interactions with the thermal cloud creates 
an effective potential for the condensate which is stiffer than the applied external potential 0]. The peak density of 
the condensate hence increases with temperature (for fixed TVq) and the critical condensate number is reduced. The 
critical number for 200 nK, as obtained from the failure to find a converged solution at larger Nq, is Nc = 1093. 
That for T — 400 nK is Nc — 1016. The temperature dependence of the critical condensate number as a function of 
temperature is shown in the inset of Fig. 2. 

It should be noted that the total number of trapped atoms, N, varies for each point in the figure. Alternatively one 
could vary T (and hence Nq) keeping A^ fixed, which would give a critical temperature for collapse. Experimentally 
evaporative cooling removes atoms. A certain total number, corresponding to the transition temperature, is reached 
at which condensation occurs. Further cooling (removal of atoms) then proceeds to a point where a second critical 
temperature (or total number) is reached, at which point the second phase transition (i.e. collapse) is observed. 
However for the experiments on ^Li the difference between the critical temperatures for BEC and collapse is extremely 
small. Cooling thus results in repeated collapse and growth, reducing A^ until a stable A'o is reached [|l3| . 

The collapse occurs because as one increases A'o, the peak noncondensate density increases due to the interactions. 
This in turn creates a tighter and tighter effective potential for the condensate, which eventually results in collapse. 
Fig. 3 shows the noncondensate density at 100 nK for a range of Nq. The dotted curve is for Aq = 50, the dashcd- 
dot-dot curve that for A^g = 1000. Note that for a large change in A'o the peak noncondensate density has changed 
relatively little. The next three curves are for A'o = 1100, 1130, and 1146, the final figure being the critical number. 
The peak noncondensate density increases rapidly over this range. This is a cooperative effect; the tighter effective 
potential reduces the frequency of the lowest collective mode. This leads to a growth in the population of this low lying 
mode, which has a density localised near the centre of the condensate, increasing the peak noncondensate density. 

The variation of the critical number with temperature is shown in the inset to Fig. 2. The rate of decrease of Nc 
with T is significantly greater in the HFB-Popov treatment than in the Hartree-Fock treatment (see Fig. 8 of Ref. 5). 
This is due to the different excitation spectra calculated in the two formalisms. In our treatment we calculate (and 
populate) the collective excitations, which include the low lying I = mode. Near collapse this mode has a much 
lower frequency than the lowest single particle excitation of the Hartree-Fock spectrum. The population of excited 
states is therefore underestimated in the Hartree-Fock treatment, and as a result, the thermal population of the state 
is lower than it is with the HFB-Popov spectrum. The noncondensate population, and hence peak density, increases 
more rapidly as a function of temperature in our calculation (c.f. inset to Fig. 3 and Fig. 7 of Ref. 5). This is what 
gives rise to the more rapid reduction in the critical number as the temperature is increased. 

In conclusion, we have presented the first self-consistent HFB-Popov calculations for a dilute gas of atoms with 
attractive effective interactions. We have studied the collective mode frequencies of such a gas and using these 
frequencies, investigated the phase transition from metastable Bose-Einstein condensate to a collapsed dense phase. 
The results from these calculations are in general agreement with previous Hartree-Fock results, but we feel that the 
HFB-Popov approach is the more appropriate one to use if only a BEC transition is assumed to take place. We find 
a significantly greater dependence of the critical number upon temperature in the HFB-Popov treatment. 

If one includes the possibility of a BCS-like pairing transition then this is not the appropriate approach as the omitted 
pair correlations (the so called anomalous average) are very important. Indeed the pair correlation term, {'ip{r)'ip(r)), 
becomes the order parameter for the BCS-like transition. The possibility of such a transition, or the existence of 
mixed phases containing both BEC and BCS macroscopic quantum states is currently under investigation. 
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FIG. 1. Noncondensate number as a function of the numerical iteration number for a trap corresponding to the Rice experi- 
ment. For A'l) = 1243 and A'o = 1244, where A'o is the number of atoms in the condensate, a converged solution is obtained with 
a self-consistent value for iV which varies by less than one part in 10* between iterations. For A'o — 1245 no stable solution can 
be found. The critical number is identified as Nc = 1244 which is the condensate number for which the I — mode frequency 
goes to zero. 

FIG. 2. Low lying mode frequencies as a function of condensate number for T = (solid), T = 200 nK (dashed), and T = 400 
nK (dotted) for a trap with a confining frequency of 150 Hz. Note the softening of the I — mode for large A*'o and the decrease 
in the critical number at which the mode goes soft as the temperature is increased. The variation of the critical number as a 
function of temperature is shown in the inset. 

FIG. 3. Noncondensate density for a range of A'o at a temperature of 100 nK in a spherical harmonic trap of frequency 150 
Hz. The curves, in increasing order of peak density, are for A'o ~ 50, 1000, 1100, 1130, and 1146. The final figure is the critical 
number. The inset shows the peak noncondensate density as a function of temperature at the critical number. 
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